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Abstract 

We introduce a model for Ce^^ impurities in cubic metals which exhibits 
competition between the Fermi-liquid fixed point of the single channel Kondo 
model and the non-Fermi-liquid fixed point of the two- and three-channel 
Kondo models. Using the non-crossing approximation and scaling theory, we 
find: (i) A possible three-channel Kondo effect between the one- and two- 
channel regimes in parameter space, (ii) The sign of the thermopower is a 
fixed point diagnostic, (iii) Our results will likely survive the introduction of 
additional and conduction states. We apply this model to interpret the 
non- Fermi liquid alloy Lai_2,.Ce2:Cu2.2Si2- 
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The common paradigm to understand the electronic properties of metals is the Landau 
Fermi-liquid picture, in which the low energy excitation spectrum has a 1:1 map to that of 
the non- interacting electron gas and concommitant constant specific heat coefficient Cei/T 
and susceptibility x{T) at low temperatures. Recent discoveries in heavy electron materials 
of non- Fermi liquid behavior in Cei/T and other properties together with non- Fermi 
liquid resistance curves in small quenched copper point contacts have stimulated renewed 
interest in the multichannel Kondo model introduced by Nozieres and Blandin P]. In this 
model an impurity with internal degrees of freedom (spin Si) couples antiferromagnetically 
to M degenerate conduction channels or bands. For M > 25/ the local spin is "overcompen- 
sated" yielding a critical ground state with a non- Fermi liquid excitation spectrum (e.g., for 
M = 2, Sj = 1/2, Cei/T ~ — ln(T) for T — |9|,p!0[|). Non-magnetic impurities with internal 



degrees of freedom generically give rise to the M = 2,Si = 1/2 model [TT||T^, and these 



models have been argued to describe uranium based heavy fermion alloy systems 
and point contact data [0. Surprisingly, to date no realistic M = 2, 5/ = 1/2 magnetic 
impurity candidates have been put forward. 

In this paper, we propose and study in detail an Anderson impurity model for a Ce^"*" 
ion in a cubic metal host which admits low temperature physics of either the M = 1, 2, or 3, 
Si = 1/2 magnetic Kondo impurity model. The model retains the lowest crystalline electric 
field (CEF) states of the and configurations. In addition to employing perturbative 
scaling arguments, we utilize the non-crossing approximation (NCA), which we calibrate by 
comparing our calculated magnetic susceptibility x(r) curves for M = 2, 3 with the exact 
results 0. We show that a large, negative thermoelectric power S{T) is a necessary condition 
for observiong the M = 2, Si = 1/2 magnetic Kondo effect. Our results will likely survive 
a more realistic treatment including all f^,P and conduction electron states. We compare 
our results to Lai_^Ce^Cu2.2Si2 @] where C(T)/T = 7(T) ~ — ln(T) and x(T) ~ — ln(T) 
(per Ce ion) are observed along with large, negative S{T). 

The original M = 2, Si = 1/2 Kondo model described an impurity magnetic moment 
with orbital conduction channel labels [§]. Heavy electron alloys containing trivalent cerium 



(Ce^+) ions with a lowest lying magnetic configuration with a doublet CEF induced 
ground state are candidate M = 2, 5*/ = 1/2 systems, provided two necessary conditions are 
met 0: 

(1) The point group of the Ce^+ ion must be cubic or hexagonal. 

(2) The ever-present M = l,Sj = 1/2 coupling (induced by /° — virtual charge fluc- 
tuations) must be smaller than the M = 2 coupling (induced by — virtual charge 
fluctuations). Equivalently, the ground state weight of must exceed that of /°. Condi- 
tions (1,2) strongly restrict observability of the magnetic M = 2, Sj = 1/2 non-Fermi liquid 
fixed point. Hence, it is not surprising that many Ce'^^-based heavy fermion materials are 
well described by the Fermi liquid fixed point of the M = 1, Sj = 1/2 Kondo model. 

Our model Anderson Hamiltonian for Ce impurities at a cubic site is developed by first 
projecting to the lowest Hund's rule LSJ multiplets for three configurations: f^{J = 0) with 
energy eo = 0, f^{J = 5/2) with energy ei, and p{J = 4) with energy e2. By taking cq = 0, 
we have €2 = 2ei + Ujf, where Uff is the on-site Coulomb repulsion. Further, we restrict 
our initial discussion to the most strongly coupled orbital angular momentum 1^ = 3, spin 
Sc = 1/2, total angular momentum = 5/2 conduction electron partial wave states. The 
conduction band density of states A^(e) is taken as lorentzian with width D. 

With the CEF turned on, the configuration with only a trivial singlet level is unaf- 
fected. We assume the CEF induces a magnetic doublet (r7) ground state in the f^J = 5/2 
multiplet which acts like an effective magnetic spin Sj=l/2 with indices a =T, i- In the 
pj = 4 multiplet, we assume a lowest lying non-magnetic doublet (Fa) which has quadrupo- 
lar or shape degrees of freedom (n="+" = stretched ion, n="-" =squashed ion). This 
roughly corresponds to the 3z^ — r^(+),x^ — ?/^(— ) orbital E doublet for transition metal 
ions in cubic symmetry. The conduction = 5/2 partial wave sextet splits into a doublet 
(Tyc) and quartet (Fgc) under the CEF. The trivial /° symmetry implies that only T^c doublet 
partial waves induce p — p mixing, p-p mixing is induced only by the Fgc = Fy^Fs partial 
waves. This quartet is a tensor product of "spin" (F7) and "orbital" (F3 ~ {x^ —y"^ , 82;^ — r^)) 
degrees of freedom. There are then three symmetry distinct "channels" of "spin states" (F7) 



from the conduction sector: one from the T^c doublet, and two from the Fgc quartet. 

Conduction impurity couphng is effected through the hybridization Hamiltonian, H^yh, 
given by 

Hhyb = K)i EcLl/° X + h.c. 

ea 

+Vu E(-l)"'''/'4J/'« X fn\ + h.c. (1) 

ena 

where Cea{cl^) and Cena(4na) annihilate (create) conduction electron partial wave states of 
energy e in the Fyc and Fgc manifolds, respectively. All Clebsch-Gordon coefficients are 
lumped into the hybridization matrix elements Vqi and V12. Although Vqi and Vu are con- 
strained by the one-body nature of the hybridization potential, we allow their free variation 
to study the competition between M = 1,2,3 parameter regimes. Similarly, we will not 
constrain ei^2 by the known values (ei ~ — 2eV, 62 ~ — 2ei ~ AeV [0]). We shall critically 
discuss these assumptions later. 

This model maps to a Kondo (exchange) form at low energy scales, using the Schrieffer- 
Wolff transformation |M. The 5/ = 1/2 impurity spin from the f^,Tj is coupled to con- 



duction states exchange integrals Ji = — 2|Voip/ei for M = l,Tic conduction states, and 
J2 = 2|l^i2p/[e2 — ei] for M = 2, Fgc conduction partial waves. When Ji = J2, an M = 3 
theory results from the combination of Tjc, Fgc states. Note that the common models with 
Uff 00 assumed give no access to the M = 2,3 fixed points! Third order scaling theory 
shows that the M = 1[2](3) fixed points obtain as Ji > J2[Ji < J2]{Ji = J2), with the 



crossover from the high temperature free moment fixed point set by the Kondo scale Tq [15 



We have studied our model with the non-crossing approximation (NCA), a self-consistent 



diagrammatic method-for further details see Refs. [[T6| , p!7| . In the NCA pseudo-particle 
Green's functions are introduced for each ionic state of the /"'^'^ configurations, and the 
self-energy equations for these propagators are self-consistently solved to second order in 
Vol, V12. The NCA yields the correct critical behavior for the over-screened SU {N)^SU{M) 
multi-channel Kondo model, with the impurity spin degeneracy |T^. The NCA produces 



pathological singularities for M = 1, but only below a scale Tp << Tq [|1^]. In particular. 



x(0),7(0) ~ 1/To as expected. Hence, the NCA is a reliable method for detailed studies of 
the competition between M = 1, 2, 3 fixed points. 

At T = 0, the NCA can be reduced to coupled non-linear differential equations for the 
inverse configuration propagators as a function of frequency which can be solved analytically 
close to a threshold energy (Eo ^ ti + 0{VQi,Vi2f) [0,|T3]. We introduce the quantity 
7c = vr[(e2 — -E'o)/ri2 + -£^0/1^01]= (7re2/ri2) which also defines a renormalized p energy €2, 
where Vij = 7tN{0)V^j. Clearly, 7^ ^ (iV(0)J2)~^ - {N{0)Ji)-\ For 7c ~ €2 > 0[< 0](= 0), 
we obtain M = 1[M = 2](M = 3) low temperature physics, in agreement with the third 
order scaling in Ji, J2 mentioned above. 

Parameters for our NCA work are listed in Table 1. For simplicity, we take F = F12 = Fqi, 
so €2 alone regulates the T — > physics. In addition, we can analytically estimate Tq (here 
ei2 = ei - £2): 

M = 1,62 > 0: keTo = e2(F/7rD)V2 exp(7rei/2F); 

M = 2,62 < 0: kBTo = (F/7r)(|e2|/Z^) ^2 exp(7rei2/2F); 

M = 3, eg = 0: fc^To = D{T/7vD)'^/^ exp(7rei/2F). 

In deriving these Tq expressions, we assume F, |ei^2| << D. For |7c| << 1, the scaling 
trajectories flow close to the M = 3 fixed point, crossing over to the M = 1,2 fixed points 
below a scale ~ |e2|^/2 H. 

Calculated x{T) curves allow us to test the NCA and develop a model phase diagram 
as shown in Fig. 1. x{T) diverges as -ln(r/ro)[(r/To)-i/^] for the M = 2[M = 3]-channel 
Si = 1/2 model as T — > [p|,p!0|, while for M = 1, x(0) ~ I/Tq at zero temperature. Our 
NCA results show approximate scaling behavior for all the M = 1, 2, 3 parameter sets, and 
agree excellently with exact Bethe ansatz results for M = 2, 3 0] as anticipated from T = 
NCA analysis [|I^ (we also find good agreement with exact results for our M = 2, 3 Cei/T 
and entropy curves [|15[). The M = 1 results are less satisfying, though x(0) ~ I/Tq in the 
NCA and negative curvature in x{T) is evident in the figure. Calculation of x"{^yT) for 
M = 2, 3 leads, as expected ||T^ to low T marginal Fermi liquid behavior 



With the NCA the electrical resistivity p(T) and thermopower S{T) may be computed 



using the Kubo formula |T6| in the dilute impurity limit, where inter- impurity correlations 



can be neglected. Assuming dominant scattering in the / = 3 partial wave sector |TH], the 



transport coefficients are determined by the integrals = / (ie(— (9//9e)r(e)e" where 
r(e) = / dkT{k,e)/4:7c is the angular averaged transport lifetime determined from the one- 
particle T-matrix defined in terms of the interconfiguration excitation propagator 1TB-|T7 . 



We obtain the expected behavior p{T)/p{0) ~ 1 - A{T/Toy/'^{l - AiT/T^f/^) power laws 
for M = 2(3) IlOlO for T < O.OSTq [Hill. 

New physics is evident in S{T) = —Ii/eTIo, with e the electron charge. Given Tq << D, 
the sign of S{T) is determined by the degree of particle-hole asymmetry in l/r(e). Figure 
2 displays the differing S(T) behaviour for the M = 1, 2, 3 regimes of our model. S(T) is 
negative definite in the M = 2, 3 regimes of the model [|1^. For M = 1, dominant /°-/^ 
mixing yields a Kondo resonance peak in 1/r just above the Fermi level, giving S(T) > 0. On 
the other hand, S{T) < arises when scattering between the and sectors dominates, 
since the overall resonance spectral weight in 1/r is shifted below the Fermi energy, giving 
stronger hole than electron scattering for T ^ 0. Hence, S{T) < for low temperatures 
when Condition (2) for the M = 2, Sj = 1/2 magnetic Kondo effect is satisfied. 

We now address the effect of omitted p and conduction states: 

(1) Contributions of excited p, Fs states We estimate J2/ Ji ~ 0.2 for realistic ri2, Fqi, £2, ei 
values and inclusion of the lowest pT^ state. However, we can make J2 > Ji within a 
spherical hybridization model (and Lorentzian density of states) by: (a) including all nine 
/^Fa states, which effectively enhances by a factor of five, and (b) reducing the conduction 
occupancy from 1.0 to 0.2 per site, qualitatively consistent with electronic structure inputs 
to analysis of spectroscopic data for CeCu2Si2(see Ref. [^(a)). 

(2) New contributions from other excited levels. Virtual charge fluctuations to excited 
triplet states suppress the one-channel coupling and induce a new channel symmetry breaking 
exchange interaction of the form H = JJ2x Sjir^Sg) where Sj, are pseudo-spin matrices 
for the F7 impurity and Fg conduction states, and the matrices from the Fg channel space 
have the forms = (-r^^) - v^r(i))/2,rS^ = (-r^^) + y3r(i))/2, and = r^^) following 



Pauli matrix conventions. Third-order scaling and strong coupling analysis show that the 
M = 2 fixed point for J2 > Ji is stable provided 2J2 > \J\, which is possible for realistic 
values of Vqi, V12, ei, £2- 

(3) Additional partial wave states. An added jc = 7/2Tsc will mix with the jc = 5/2, Fgc 
quartet resulting in new exchange terms. A unitary transformation to bonding and non- 
bonding quartets yields two identical, decoupled M = 2, Si = 1/2 equations with enhanced 
coupling strength in the bonding channel. The non-bonding coupling has a zero bare value, 
and so is practically irrelevant. 

Turning to the alloy, Lai_2;Ce^Cu2.2Si2, it is found that x(T),7(T) ~ — ln(T) for x = 
0.1 1^. We thus examine whether it meets the necessary Conditions (1,2) for the M = 
2, Sj = 1/2 fixed point. First, the Ce ions sit in a pseudo-cubic environment with a doublet 
ground state well known from neutron scattering [|19| . The best superconducting samples of 



CeCu2.2Si2 have isotropic low temperature susceptibilities |21[|. Hence, condition (1) appears 



to be satisfied. Using data from Ref. |^ and the cubic effective moment to estimate the 
Landau- Wilson ratio R given by 

'7T'kUdx{T)/dlnT) 



R 



(2) 



filjf{dj{T)/d\nT)\ 

we obtain R =2.7(1), in excellent agreement with the theoretical M = 2, Sj = 1/2 value 
R = 8/3 [|,|l3]- The uncertainty reflects measurement accuracy from the plots of Ref. 0. 
Further support arises from the enhancement of 7(T) in applied magnetic fleld which is 
qualitatively (though not quantitatively) consistent with the M = 2, Sj = 1/2 Bethe-Ansatz 
results 0; in contrast, for the M = 1, Si = 1/2 model, 7(T) drops with applied fleld. The 
resistivity p(T) does not display clean T^/^ behavior, but has been measured down to only 
T ^ To/10=1.2K IQ; our work shows that the T^/^ behavior sets in below about O.OSTq 



15,18 



Turning to condition (2), we consider the x = 1 S{T) data for which a sign change 
occurs at 70K, well above To ~ lOK, eventually reaching a large negative value of about 
—20 — 30yuV/K, which is in reasonable agreement with our M = 2 curves. We observe that: 



(i) For X = the thermopower is positive, positively identifying the sign change with the 
Ce ions, (ii) Coherent lattice quasiparticle effects can only be expected for T << Tg ~ lOi^' 
(in other Ce based heavy fermion materials, sign changes typically occur well below Tq 
with magnitudes of only a few /iV/K). (iii) The sign change occurs well below the CEF 
excitation energy of 360K pT| , P^ , so is not due to CEF effects (although the positive S(T) 
values for higher T are likely due to scattering off the excited CEF quartet, which generates 
a Kondo resonance above the Fermi level), (iv) The magnitude is easily explained only 
with a resonant scattering mechanism near the Fermi level. Observations (i)-(iv) strongly 
suggest that the negative thermopower arises from the single ion Kondo effect with — 
fluctuations dominating the virtual exchange processes, hence satisfying necessary condition 



(2), which is further supported by analysis of high energy spectroscopies pO|. Negative 
S{T) are measured at T = 20K for many Ce based 1-2-2 structures with high unit cell 
volumes, while positive S{T) at T = 20K are seen for low volumes |2^. Indeed, increasing 
pressure can tune from M = 2 through M = 3 to M = 1 fixed points since pressure initially 
destabilizes p relative to /° |T^. 

In summary, we believe the Lai_^Ce^Cu2.2Si2 system is a promising M = 2, 5*/ = 1/2 (or 
M = 3 under pressure) Kondo candidate which should be studied further, particularly be- 
cause magnetic Kondo impurities are more readily studied (they couple directly to magnetic 
field and neutrons). 
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TABLES 



Set 


M 


r/D 


ei/D 


e2/D 


kBT^/D 


1 


2 


0.2 


-0.4 


-0.1 


.19081E-02 


2 


2 


0.2 


-0.37 


-0.07 


.15964E-02 


3 


2 


0.2 


-0.35 


-0.05 


.13492E-02 


4 


3 


0.2 


-0.3 


0.0 


.15224E-02 


5 


3 


0.2 


-0.4 


0.0 


.69413E-03 


6 


1 


0.2 


-0.3 


0.05 


0.11957E-02 


7 


1 


0.2 


-0.3 


0.07 


0.16740E-02 


8 


1 


0.2 


-0.3 


0.10 


0.23914E-02 



TABLE I. Model parameters for the Ce impurity. This set of model parameters covers the 
single, three, and two channel Kondo regimes. F = Fqi = is the hybridization strength for 
both — and — mixing, respectively. 



Figure 1. Scaling behavior for x{T) of model Ce'^+ impurity. Results are shown 
for parameter regimes of the model (Eq. (1)) giving M = 1,2,3 fixed points (see Table 1). 
Points: NCA calculations; Solid Lines: Exact Bethe-Ansatz results [^J (with Tq = Tk /0.3 
for M = 2, Tk from Ref. 0). For convenience, we multiply x{T) by 2.0 for M = 2. The 
ground state phase diagram for the model in exchange coupling constant parameter space 
is drawn in the inset, where gi = N{0)Ji, N{0) being the conduction band density of states 
at the Fermi energy. The solid diagonal line is for M = 3. 

Figure 2. Thermopower S(T). The thermopower is positive at low T for M = 1, 
strongly negative for M = 2, and weakly negative for M = 3. Dominant /° — /^(/^ — /^) 
virtual charge fluctuations give positive (negative) S{T). Referring to the Table 1: the dash- 
dotted line is calculated with parameter set 8, the dashed-dot-dotted line is calculated with 
parameter set 1, and the solid line is calculated with parameter set 4. 
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